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Abstract: Let {F„ : n ^ 1} be a normalized sequence of random variables in some fixed Wiener 
chaos associated with a general Gaussian field, and assume that — >• E[N^] = 3, where 
is a standard Gaussian random variable. Our main result is the following general bound: there 
exist two finite constants c, C > such that, for n sufficiently large, c x max(|£'[F^]|, E[F^] — 3) ^ 
d{Fn,N) xmsix{\E[F^]\,E[F^]-3), where d{Fn,N) = sup \E[h{Fn)] - E[h{N)]\, and h runs 
over the class of all real functions with a second derivative bounded by 1. This shows that the 
deterministic sequence max(|£^[F^]|, — 3), n ^ 1, completely characterizes the rate of con- 
vergence (with respect to smooth distances) in CLTs involving chaotic random variables. These 
results are used to determine optimal rates of convergence in the Breuer-Major central limit the- 
orem, with specific emphasis on fractional Gaussian noise. 
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1 Introduction and main results 

Let X = {X{h) : h £ i^} be an isonormal Gaussian process (defined on an adequate space 
(ri, J^, P)) over some real separable Hilbert space Sj, fix an integer q ^ 2, and let {F„ : n ^ 1} be 
a sequence of random variables belonging to the g'th Wiener chaos of X (see Section [2. II for precise 
definitions). Assume that = 1 for every n. In recent years, many efforts have been devoted 

to the characterization of those chaotic sequences {Fn} verifying a Central Limit Theorem (CLT), 
that is, such that Fn converges in distribution to ~ ^(0,1) (as n — )■ oo), where ^(0,1) 
denotes a centered Gaussian law with unit variance. An exhaustive solution to this problem was 
first given by Nualart and Peccati in [17) . in the form of the following "fourth moment theorem". 
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Theorem 1.1 (Fourth Moment Theorem — see |17j ) Fix an integer q ^ 2, and consider a 
sequence of random variables {F„ : n ^ 1} belonging to the qth Wiener chaos of X and such that 
E[F^] = 1 for all n ^ 1. Then, as n ^ oo, Fn converges in distribution to N ^ =^(0, 1) if and 
onlyifE[F^\^E[N^\ = 'i . 

Note that Theorem 11.11 represents a drastic simplification of the usual method of moments 
and cumulants, as described e.g. in |18| . Combining the so-called Stein's method for normal 
approximations (see |H [13] , as well as Section 13.11 below) with Malliavin calculus (see [51 [16] , as 
well as Section 12. 2p , one can also prove the forthcoming Theorem 11.21 providing explicit upper 
bounds in the total variation distance. We recall that the total variation distance dTv{F-,G) 
between the laws of two real- valued random variables F, G is defined as 



dTv{F,G) 



sup 



\P[F G A]-P[G G ^]| 



where the supremum runs over the class of all Borel sets ^ C M. Given a smooth functional of 
the isonormal process X, we shall also write DF to indicate the Malliavin derivative of F (thus 
DF is a i^-valued random element - see again Section [2.21 for details). 

Theorem 1.2 (Fourth Moment Bounds — see [9|, I15| ) Fix q ^ 2, let F be an element of the 
qth Wiener chaos of X with unit variance, and let N ~ ^(0,1). The following bounds are in 
order: 



dTv{F,N) ^ 2 



E 



^--\\DF\\l 



^2W^(i^[F4]-3). 



3<Z 



(1.1) 



Remark 1.3 1. The two inequalities in (jl.ip were discovered, respectively, in [9] and |15) . 
Using the properties of the Malliavin derivative DF (see Section 12.21 below) , one sees im- 
mediately that 



E 



1 - -\\DF\ 



Var ( -||L>F|| 



2. One can prove the following refinement of the second inequality in (jl.ip (see |1H Lemma 
3.5]): for every random variable F belonging to the gth Wiener chaos of X and with unit 
variance 



1, 



q-l 



Var ( - \\DF\\% ) ^ ^-^{E[F^] - 3) ^ (g - l)Var ( -\\DF 



1, 



(1.2) 



3. Theorem II. 21 implies that, not only the condition -^[-F^] — )• 3 is necessary and sufficient for 
convergence to Gaussian, as stated in Theorem ll.il but also that the sequences 



^(n) := ^E[F^] - 3 and 7( 



n 



'Var ( -^WDFnWl 



n > 1, 



(1.3) 



bound from above (up to a constant) the speed of convergence of the law of Fn to that of 
N in the topology induced by dTV- 
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4. If one replaces the total variation distance with the Kolmogorov distance or with the Wasser- 
stein distance (see e.g. [9l [13] for definitions), then the bounds (jl.ip hold without the 
multiplicative factor 2 before the square roots. 

5. When E[F] = and E[F'^] = 1, the quantity — 3 coincides with the fourth cumulant 
of F, see Definition 13.31 One can also prove that, if F belongs to a fixed Wiener chaos and 
has unit variance, then > 3 (see |17|). 

6. Throughout the paper, in order to simplify the notation, we only consider sequences of ran- 
dom variables having unit variance. The extension to arbitrary sequences whose variances 
converge to a constant can be deduced by a straightforward adaptation of our arguments. 

A natural problem is now the following. 

Problem 1.4 Assume that {Fn} is a unit variance sequence belonging to the qth Wiener chaos of 
the isonormal Gaussian process X. Suppose that Fn converges in distribution to N ^ ^(0, 1) and 
fix a distance do{-, ■) between the laws of real-valued random variables. Can one find an explicit 
optimal rate of convergence associated with the distance do ? 

The notion of optimality adopted throughout the paper is contained in the next definition. 

Definition 1.5 Assume that, as n — t- go, Fn converges in distribution to N, and fix a generic 
distance do{-, ■) between the laws of real- valued random variables. A deterministic sequence {(p{n) : 
n ^ 1} such that <p{n) ^ is said to provide an optimal rate of convergence with respect to do if 
there exist constants < c < C < oo (not depending on n) such that, for n large enough, 

c<;M^<a (1,4) 

The problem of finding optimal rates is indeed quite subtle. A partial solution to Problem 11.41 
is given by Nourdin and Peccati in |10| . In this reference, a set of sufficient conditions are derived, 
ensuring that the sequences /3(n),7(n) in (jl.3p yield optimal rates for the distance do = dxv- In 
particular, these conditions involve the joint convergence of the two-dimensional vectors 

l-q-^\\DFn\\l \ ^ ^ 

Fn, r-\ ^ ' n ^ 1. (1.5) 



7(n) 

The following statement constitutes one of the main finding of |10) (note that the reference [10 



only deals with the Kolmogorov distance but, as far as lower bounds are concerned, it is no more 
difficult to work directly with dTv)- 

Tiieorem 1.6 (See [lOj) Let {Fn} be a unit variance sequence belonging to the qth Wiener chaos 
of X, and suppose that, as n ^ oo, Fn converges in distribution to N ^ ^(0,1). Assume 
moreover that the sequence of two-dimensional random vectors in M.5\) converges in distribution 
to a Gaussian vector {Ni,N2) such that E[N^] = E[N^] = 1, and E[NiN2] =: p e (-1, 1). Then, 
for every z € M.' 

/2 

7(?i)^^[P(F„ ^ z) - P(iV < z)] ^ ^(z^ - 1)1^^ asn^oo. (1.6) 

3 \/2tt 



In particular, if p ^ the sequences f3{n),j{n) defined in U.3\) provide optimal rates of convergence 
with respect to the total variation distance dTVi in- the sense of Definition \l.[A 
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As shown in |101 Theorem 3.1] the conditions stated in Theorem 11.61 can be generahzed to 
arbitrary sequences of smooth random variables (not necessarily belonging to a fixed Wiener 
chaos). Moreover, the content of Theorem 11.61 can be restated in terms of contractions (see \10\ 
Theorem 3.6]) or, for elements of the second Wiener chaos of X, in terms of cumulants (see |10^ 
Proposition 3.8]). 

One should note that the techniques developed in [10] also imply analogous results for the 
Kolmogorov and the Wasserstein distances, that we shall denote respectively by dpcol and dw- 
However, although quite flexible and far-reaching, the findings of |10j do not allow to deduce a 
complete solution (that is, a solution valid for arbitrary sequences {Fn} in a fixed Wiener chaos) 
of Problem 11.41 for either one of the distances dxvj dxol and dyy. For instance, the results of |10j 
provide optimal rates in the Breuer-Major CLT only when the involved subordinating function 
has an even Hermite rank, whereas the general case remained an open problem till now - see \10\ 
Section 6]. 

The aim of this paper is to provide an exhaustive solution to Problem 11.41 in the case of 
a suitable smooth distance between laws of real-valued random variables. The distance we are 
interested in is denoted by d{-,-), and involves test functions that are twice differentiable. The 
formal definition of d is given below. 

Definition 1.7 Given two random variables F,G with finite second moments we write d{F,G) 
in order to indicate the quantity 

d{F, G) = sup \E[h{F)] - E[h{G)] I , 

where U stands for the set of functions /i : M — > M which are (that is, twice differentiable and 
with continuous derivatives) and such that ||/i"||oo ^ 1- 

Observe that d{-,-) defines an actual distance on the class of the distributions of random 
variables having a finite second moment. Also, the topology induced by d on this class is stronger 
than the topology of the convergence in distribution, that is: if d{Fn,G) — t- 0, then Fn converges 
in distribution to G. 

The forthcoming Theorem ll.91 and Theorem II. Ill contain the principal upper and lower bounds 
proved in this work: once merged, they show that the sequence 

max{\E[F^]lE[F^]-3}, n ^ 1, (1.7) 

always provides optimal rates for the distance d, whenever {Fn} lives inside a fixed Wiener chaos. 
As anticipated, this yields an exhaustive solution to Problem 11.41 in the case do = d. One should 
also note that the speed of convergence to zero of the quantity (|1.7|) can be given by either one 
of the two sequences {|£'[-F^]|} and — 3}; see indeed CoroUarv 16.81 for explicit examples of 

both situations. 

Remark 1.8 Let {Fn : n ^ 1} be a sequence of random variables living inside a finite sums of 
Wiener chaoses. Assume that — ?• 1 and that Fn converges in distribution to ~ ^(0, 1). 
Then, the hypercontractivity property of the Wiener chaos (see e.g. [T] Chapter V]) imply that 
E[Fn] — )■ E[N^] for every integer A; ^ 3. In particular, one has necessarily that E[Fn] — )• 0. 
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Theorem 1.9 (Upper bounds) Let N ~ ^(0, 1) be a standard Gaussian random variable. 
Then, there exists C > such that, for all integer q ^ 2 and all element F of the qth Wiener 
chaos with unit variance, 

d{F,N) ^ C max {\E[F^,E[F'^] - 3}. (1.8) 

Remark 1.10 1. In the statement of Theorem II .91 we may assume without loss of generahty 
that N is stochastically independent of F. Then, by suitably using integration by parts and 
then Cauchy-Schwarz (see e.g. |14| Theorem 3.2]), one can show that, for every h £U (see 



Definition [LT]) . 

\E[h{F)] - E[h{N)]\ 



E h"{Vl^tF + ViN)(^-^\\DF\\l^ 



dt 



Since (jl.2p is in order, one sees that the inequality (jl.9p does not allow to obtain a better 
bound than 



d{F, N) C^E[F^] - 3, 
which is not sharp in general, compare indeed with (jl.Sp . One should observe that the 



rate ^JE\F^\^^ may happen to be optimal in some instances, precisely when E[F^\ and 
^JE\F^\^^ have the same order. In the already quoted paper |10) one can find several 
explicit examples where this phenomenon takes place. 

2. Let F be an element of the qth Wiener chaos of some isonormal Gaussian process, and 
assume that F has variance 1. It is shown in [111 Proposition 3.14] that there exists a 
constant C, depending only on q, such that [^^[F^]! ^ C-\/E[F'^\ — 3. Using this fact, one 
can therefore obtain yet another proof of Theorem 1 1 . 1 1 based on the upper bound (|1.8p . 



Theorem 1.11 (Lower bounds) Fix an integer 2 and consider a sequence of random vari- 
ables {Fn : n ^ 1} belonging to the qth Wiener chaos of some isonormal Gaussian process and 
such that E[F^] = 1. A ssume that, as n ^ oo, Fn converges in distribution to N ^ ^(0,1). 
Then there exists c > (depending on the sequence {Fn}, but not on n) such that 

d(F„,iV) ^ c X max{|^[F^]|,^[F^] - 3), n^l. (1.10) 

Our proofs revolve around several new estimates (detailed in Section Uj), that are in turn based 
on the analytic characterization of cumulants given in |12] . Also, a fundamental role is played by 
the Edgeworth-type expansions introduced by Barbour in [T]. 



1.1 Plan 

The paper is organized as follows. Section [2] is devoted to some preliminary results of Gaussian 
analysis and Malliavin calculus. Section |3] deals with Stein's method, cumulants and Edgeworth- 
type expansions. Section 3] contains the main technical estimates of the paper. Section [5] focuses 
on the proofs of our main findings, whereas in Section |6] one can find several applications to the 
computation of optimal rates in the Breuer-Major CLT. 
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2 Elements of Gaussian analysis and Malliavin calculus 



This section contains the essential elements of Gaussian analysis and Malliavin calculus that are 
used in this paper. See the classical references [8l [16] for further details. 

2.1 Isonormal processes and multiple integrals 

Let Sj he a. real separable Hilbert space. For any q ^ 1, we write i^®"? and i^®'^ to indicate, 
respectively, the gth tensor power and the qth symmetric tensor power of i^; we also set by 
convention S)®^ = i^®'^ = M. When = L'^(A,A,fi) =: L'^{fi), where /i is a a-finite and non- 
atomic measure on the measurable space [A, A), then i^®"^ = L'^ [A'^ , A'^ , fi'^) =: and 
i^®? = Ll{A'' , A'^ , fi'') := LKi^i''), where L'^ifi'') stands for the subspace of L?{^'^) composed of 
those functions that are /i'^-almost everywhere symmetric. We denote by X = {X{h) : h E 
i^} an isonormal Gaussian process over i^. This means that X is a centered Gaussian family, 
defined on some probability space {Q,J^,P), with a covariance structure given by the relation 
E [X{h)X{g)] = {h,g)sj- We also assume that = cr{X), that is, is generated by X. 

For every q ^ 1, the symbol Tig stands for the qth Wiener chaos of X, defined as the closed 
linear subspace of L'^{Q,,J^,P) =: L^(Q) generated by the family {Hq{X{h)) : h £ Sj, = 1}, 
where Hq is the g'th Hermite polynomial given by 

Hq{x) = {-iye^ — {e-^). (2.11) 

We write by convention T-Lq = M. For any q ^ 1, the mapping Iq[h®'^) = Hq{X[h)) can be 
extended to a linear isometry between the symmetric tensor product S^®"^ (equipped with the 
modified norm IHIj-,®?) and the gth Wiener chaos Hq. For q = 0, we write Io{c) = c, c G M. 
A crucial fact is that, when fj = L^(/x), for every / S i^®'^ = L^(/x'^) the random variable Iq{f) 
coincides with the q-fold multiple Wiener-Ito stochastic integral of / with respect to the centered 
Gaussian measure (with control /i) canonically generated by X (see |16l Section 1.1.2]). 

It is well-known that L^(0) can be decomposed into the infinite orthogonal sum of the spaces 
Hq. It follows that any square-integrable random variable F G L^(J7) admits the following Wiener- 
Ito chaotic expansion 

oo 

<?=o 

where /q = and the fq € io®'', g ^ 1, are uniquely determined by F. For every g ^ 0, 

we denote by Jq the orthogonal projection operator on the qth Wiener chaos. In particular, if 
F e L^{Q) is as in ff2J2]l . then JqF = Iq{fq) for every g ^ 0. 

Let {cfc. A; ^ 1} be a complete orthonormal system in i^. Given / G io®^ and g G i^®'', for 
every r = 0, . . . ,p A the contraction of / and g of order r is the element of i5®(p+9-2r) defined 

by 

oo 

f®rg= ^ (/, (g) ... (g) ei^)^®r (g) (5,6,1 (g) ... (g) ei^)j5®r. (2.13) 

il,...,jr = l 

Notice that the definition of / (g^, gi does not depend on the particular choice of {cfc, k ^ 1}, and 
that f ®r g is not necessarily symmetric; we denote its symmetrization by f®rg G j^Q(p+'?~2r)^ 
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Moreover, f®Qg = f®g equals the tensor product of / and g while, for p = q, f ®qg = {f,g)sjSS"i- 
When S) = L'^{A,A,fi) and r = 1, ...,p A q, the contraction / Cg)^ 5 is the element of L'^{fx^^'^~'^^) 
given by 

f i^r g{xi, ■.■,Xp+g-.2r) (2.14) 
/(Xi,... )g{Xp_r+l, ■■■,Xp+g-.2r, ^1, . , ar)d^{ai) . ..dfl{ar) ■ 



It can also be shown that the following multiplication formula holds: if / G Sj®^ and g G i^®'', 
then 

Ip{f)Ig{g) = Ilj\r) (^J W2r(/§r.ff). (2.15) 



2.2 Malliavin operators 

We now introduce some basic elements of the Malliavin calculus with respect to the isonormal 
Gaussian process X. Let S be the set of all cylindrical random variables of the form 

F = g{X{cPi),...,X{4>n)), (2.16) 

where n ^ 1, g : — ?• R is an infinitely diff'erentiable function such that its partial derivatives 
have polynomial growth, and (pi £ Sj, i = 1, . . . ,n. The Malliavin derivative of F with respect to 
X is the element of L^(r2,io) defined as 

i=l 

In particular, DX{h) = h for every h £ Sj. By iteration, one can define the mth derivative D^F, 
which is an element of L^(f], i^®™), for every m ^ 2. For m ^ 1 and p ^ 1, denotes the 

closure of S with respect to the noriTi II * ||777,^r), defined, by the relation 



m 

= E[m + Y,E[\\D^F 



m 

p 

i=l 



We often use the notation := n,„^i flp^i 

Remark 2.1 Any random variable Y that is a finite linear combination of multiple Wiener-Ito 
integrals is an element of D°°. Moreover, if y 7^ 0, then the law of Y admits a density with respect 
to the Lebesgue measure - see |21j . 



The Malliavin derivative D obeys the following chain rule. If ip : — t- IR is continuously 
differentiable with bounded partial derivatives and if F = (Fi, . . . ,Fn) is a vector of elements of 
B^'^ then ip{F) E B^'^ and 

D^(P) = T.Ap^DF,. (2.17) 

1=1 * 
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Remark 2.2 By approximation, it is easily checked that equation (|2.17p continues to hold in the 
following two cases: (i) Fi G D°° and has continuous partial derivatives with at most polynomial 
growth, and (ii) Fi G D^''^ has an absolutely continuous distribution and 99 is Lipschitz continuous. 

Note also that a random variable F as in (|2.12p is in O^'"^ if and only if X]^]^ g|| Jq-F||^2(j^) < 00 
and, in this case, = '^'^i Q\\JqF\\'j^2(^Qy H = L'^ {A, A, fi) (with /x non-atomic), then 

the derivative of a random variable F as in (|2.12|) can be identified with the element of L'^{A x Q) 
given by 

00 

D,F = Y,Qlg-iifqi;x)), xeA. (2.18) 

q=l 

We denote by 6 the adjoint of the operator D, also called the divergence operator. A ran- 
dom element u G L'^{Q,S)) belongs to the domain of 6, noted Dom(5, if and only if it verifies 
\E{DF,u)fj\ ^ Cm ||-F||j;^2(Q) for any F £ D^'^, where Cu is a constant depending only on u. If 
u G Dom5, then the random variable 6{u) is defined by the duality relationship (customarily 
called integration by parts formula) 

E[F6{u)]=E[{DF,u),,], (2.19) 

which holds for every F G O^''^. 

The operator L, defined as L = X^^g ~1-^qj ^^'^ infinitesimal generator of the Ornstein- 
Uhlenbeck semigroup. The domain of L is 



DomL = {F£ L'^in) : || J^FH^a^j^^ < 00} 

q=l 



^2,2 



There is an important relation between the operators D, 5 and L. A random variable F belongs 
to D^'^ if and only if F G Dom {6D) (i.e. F G D^'^ and DF G Dom5) and, in this case, 

6DF = -LF. (2.20) 

For any F G L'^{Q), we define L~^F = Y^'^^i -^Jq{F). The operator L^^ is called the 
pseudo-inverse of L. Indeed, for any F G L^(f]), we have that L~^F G DomL = B^'^, and 

LL~^F = F - E{F). (2.21) 
The following result is used throughout the paper. 
Lemma 2.3 Suppose that H G B^'^ and G G L'^{9). Then, L'^G G O^'"^ and 

E[HG] = E[H]E[G] + E[{DH, -DL-^G)sj]. (2.22) 
Proof By (f2:20]l and <^2B), 

E[HG] - E[H]E[G] = E[H{G - E[G])] = E[H x LL^^G] = E[H6{-DL-^G)], 
and the result is obtained by using the integration by parts formula ()2.19p . |--| 
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3 Stein's equations and cumulants 



In order to prove our main results, we shall combine the integration by parts formula of Malliavin 
calculus, both with a standard version of the Stein's method for normal approximations (see [3] for 
an exhaustive presentation of this technique) and with a fine analysis of the cumulants associated 
with random variables living in a fixed chaos. One of our main tools is an Edgeworth-type 
expansion (inspired by Barbour's paper [T]) for smooth transformations of Malliavin differentiable 
random variables. These fundamental topics are presented in the three subsections to follow. 



3.1 Stein's equations and associated bounds 

Let ~ ^(0, 1) be a standard Gaussian random variable, and let /i : M — )■ M be a continuous 
function. 

Remark 3.1 In the literature about Stein's method and normal approximation, it is customary 
at this stage to assume that h is merely Borel measurable. However, this leads to some technical 
issues that are not necessary here. See e.g. |13| Chapter 3]. 

We associate with h the following Stein's equation: 

h{x) - E[h{N)] = f'{x) - xf{x), X G M. (3.23) 

It is easily checked that, if E\h{N)\ < oo, then the function 

Mx) = e^'/2 r (/^(y) _ E[h{N)])e-y'/''dy, x G M, (3.24) 

is the unique solution of ()3.23p verifying the additional asymptotic condition 
hrn h{x)e--y^=0. 

x— >±oo 

In this paper, we will actually deal with Stein's equations associated with functions h that 
are differentiable up to a certain order. The following statement (proved by Daly in [5]) is an 
important tool for our analysis. Throughout the following, given a smooth function <^ : M — )■ M, we 
shall denote by g^''\ k = 1,2, the kth derivative of g; we sometimes write g' = g^^\ g" = g^'^\ 
and so on. 



Proposition 3.2 Let the previous notation prevail, fix an integer k ^ 0, and assume that the 
function h is {k + l)'times differentiable and such that h^^^ is absolutely continuous. Then, fh is 
{k + 2) -times differentiable, and one has the estimate 

||/f+')|U^2||/i(^+i)||oo. (3.25) 
Moreover, the continuity of implies the continuity of f^'^'^^ . 

Proof. The first part, i.e., inequality (|3.25|) . is exactly Theorem 1.1 of [5], whereas the transfer of 
continuity is easily checked by induction and by using ()3.23p . |--| 
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3.2 Cumulants 

We now formally define the cumulants associated with a random variable. 

Definition 3.3 (Cumulants) Let F be a real-valued random variable such that E\F\^ < oo 
for some integer m ^ 1, and define (t)F{t) = -^[e**^], t G M, to be the characteristic function of F. 
Then, for j = 1, ...,m, the jth cumulant of F, denoted by Kj{F), is given by 



K,{F) = i-iy — logMt)\t=o. (3.26) 



Remark 3.4 The first four cumulants are the following: ki{F) = E[F], K2{F) = ElF"^]- ElF]"^ = 
Var(F), K3(F) = E[F^] - 3E[F'^]E[F] + 2E{Ff, and 

K4(F) = E[F'^] - 3E[F]E[F^] - 3E[F'^f + UElF]"^ ElF"^] - QE[Ff. 

In particular, when E[F] = one sees that K3(F) = ElF'"^] and K4(F) = ElF"^] - SElF"^]^. 

The reader is referred to |18| Chapter 3] for a self-contained presentation of the properties of 
cumulants and for several combinatorial characterizations. The following relation (proved e.g. in 
\1'2\ Proposition 2.2]) shows that moments can be recursively defined in terms of cumulants (and 
vice-versa): fix m = 1,2..., and assume that E\F\"^^^ < oo, then 

^ Y,[)^^s+i{F)E[F^~']. (3.27) 



s=0 



We now want to characterize cumulants in terms of Malliavin operators. To do so, we need 
the following recursive definition (taken from |12)). 



Definition 3.5 Let F G B°°. The sequence of random variables {Tj{F) : j ^ 0} C B°° is 
recursively defined as follows. Set ro(-F) = F and, for every j ^ 1, 

TjiF) = {DF,-DL-^Tj^,{F)),,. 

Note that each Tj{F) is a well-defined element of D°°, since F is assumed to be in B°° - see [HI 
Lemma 4.2(3)] 

For instance, one has that Ti{F) = {DF,—DL^^F)f,. The following statement provides two 
explicit relations ( ()3.28p and p.29p ) connecting the random variables Tj{F) to the cumulants of 
F. Equation (|3.28p has been proved in |12| Theorem 4.3], whereas (j3.29p is new. 

Proposition 3.6 Let F G B°°. Then F has finite moments of every order, and the following 
relation holds for every s ^ 0: 

Ks+i{F) = slE[rs(F)]. (3.28) 

// moreover E{F) = then, for every s ^ 1, 



K,+2(F) = i(s + l)!E 



(.-1)! 



(3.29) 
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Proof. In view of \12\ Theorem 4.3], we have only to prove (|3.29p . Applying Lemma 12.31 in the 
special case H = and G = Tg^iiF), and using the relation DF^ = 2FDF, one deduces that 

E[F^rs^i{F)] = E[F^]E[rs^i{F)] + 2£;[Fr,(F)]. 

Now apply Lemma 12.31 in the case H = F and G = Vs {F) : exploiting the fact that F is centered 
together with (j3.28p . we infer that 

E[FT,{F)] = E[F{T,{F) - E\T,{F)])] = E\T,^,{F)] = j^^- (3.30) 

Since ()3.28p implies that (s — l)!i?[rs_i(-F)] = Ks{F), the conclusion follows. 



Remark 3.7 L Relation ()3.28p continues to hold under weaker assumptions on the regularity 
of F. See again |12| Theorem 4.3]. 

2. Relation ()3.28p generalizes the following well-known fact: if F G D-^'^, then ri(F) = 
{DF, -DL-^F)rj is in L'^{n) and 

Var(F) = E[ri(F)]. (3.31) 

The following statement provides an explicit expression for Ts{F), s ^ 1, when F has the 
form of a multiple integral. 

Proposition 3.8 (see |12| . formula (5.25)) Let q ^ 2, and assume that F = Iq{f) with f G 
S^®'' . Then, for any s ^ 1, we have 

q [sq-2ri-...-2rs-i]/\q 

^s{F) = XI--- Cq{ri,...,rs)l{r^<q}...l^ri + ...+Ts-^<'-^} (3-32) 

ri=l rs=l 

Xl(^+l)g_2,.,_..._2r, {{■■■{f®rj)®r2f) ■ ■ ■ f)®rj), 

where the constants Cg(ri, . . . , rs-2) o-re recursively defined as follows: 

2 



Cq{r) = q{r - 1)! 
and, for a ^ 2, 



q-1 
r — 1 



aq-2ri- ... - 2ra-i - 1\ / 9 - 1 



Cq{ri,...,ra)=q{ra-l)\[ j (^^^ _ ^ )cg(ri, . . . ,r„_i). 



Remark 3.9 By combining (|3.28p with ()3.32p . we immediately get a representation of cumulants 
that is alternative to the one based on 'diagram formulae'. See |121 Theorem 5.1] for details on 
this point. 
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3.3 Assessing Edgeworth-type expansions 

The following Edgeworth-type expansion also plays a crucial role in the following. 

Proposition 3.10 Let F be an element o/B°°. Then, for every M ^ 1 and every function 
/ : M — 7- M that is M times continuously differentiahle with derivatives having at most polynomial 
growth, we have 

E[Ff{F)] = ^^^^E[f(^\F)] + E[rM{F)f^''HF)]. (3.33) 
.=0 ^■ 

Proof. Using twice Lemma 12. 3| first in the case H = F and G = f{F) and then in the case 
F = Ti{F) and G = f'{F), we deduce that 

E[Ff{F)] = E[F]E[f{F)]+E[f'{F)Ti{F)] 

= E[F]E[f{F)] + E[f'{F)]E[r,iF)] + E[f {F)r2{F)], 

where we have used the chain rule (j2.17p as well as Remark 12.21 Therefore, (j3.33p holds for 
M = 1,2 (see also (|3.31|) ). The case of a general M follows immediately from an induction 
argument and by using ()3.28p . ^ 

The following statements contain two important consequences of (|3.33p . They will be used in 
order to prove our main findings. 

Corollary 3.11 Let N ~ ^(0, 1) and fix F € B°° such that E[F] = 0, ElF"^] = K2{F) = 1. For 
M ^ 2 , let h : ^ R be (M — 1) times continuously differentiable with bounded derivatives, and 
define fh according to \3.24^ . Then, 

^2\\h'^^'-^^UE\TM{F)\. (3.34) 

Proof. From Proposition l3.2| we deduce that the function f^ is M-times continuously differentiable 
and that, for k = 2,...,M, \\f]^''\\oc ^ 2||/i(*^-i) ||oo. Using a Taylor expansion, we deduce that 
has at most polynomial growth. It follows that (|3.33p can be applied to the function fh, and the 
conclusion is obtained from the relation E[h{N)] - E[h{F)] = E[Ffh{F)] - E[f'^{F)]. 



E[h{N)] - E[h{F)] 



E '-^E[fi:\F)] 



s\ 



Corollary 3.12 Let N ~ ^(0, 1) and fix F £lD)°^ such that E[F] = and E[F'^] = K2{F) = 1. 
Let /i : M — 7- M &e twice continuously differentiable and such that ||/i"||oo ^ 1> (^nd define fh 
according to \3.24^ . Then, 

\E[h{F)] - E[h{N)] \ i^K\E[F'^]\+2E\T^{F)\- (3.35) 

where K := 1 + E[\F\]. 
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Proof. We first observe that E[h{F)]- E[h{N)] = E[h{F)]- E[h{N)], where h{x) = h{x)-h{0)- 
h'{0)x, so that we can assume without loss of generahty that h{0) = h'{0) = 0. Thus, because 
\\h"\\oo ^ 1, we get that \h{x)\ ^ ^ and \h'{x)\ ^ for ah x G M, while \E[h{N)]\ ^ i. It follows 
from ([3:2^11 that 

Next, Proposition 13 . 21 shows that is thrice continuously differentiable with ||//^"||oo ^ 2 ||/i"||oo ^ 
2. On the other hand, for all j; G R, 

f'^,{x) = xfh{x) + h{x)-E[h{N)], 
f'llix) = fi,{x) + xfUx) + h'{x). 

Consequently, /^'(O) = fh{0) and 

\f'ax)\ ^ |A(0)| + \f'^{x) - fm ^ 2 + ll/niookl ^ 2 + 2|x|. (3.36) 

We deduce that \E[fll{F)]\ ^ 2K. Applying ([3:33]) to fh in the case M = 3 yields therefore 

E[FMF)] = EifUF)] + l^E[mF)]E[F^] + E[f'^' {F)r,{F)], (3.37) 
implying in turn that 

\E[hiF)] - E[h{N)]\ ^ ^\EmF)]\\E[F^]\ + |/riooi^|r3(F)|, 
from which the desired conclusion follows. 

□ 



Remark 3.13 1. The idea of bounding quantities of the type 

M-l , . 

E[Ff{F)]-Y,^^±^E[f(^HF)] , 

8 = 

in order to estimate the distance between F and N ~ ^{0,1), dates back to Barbour's 
seminal paper [1]. Due to the fact that F is a smooth functional of a Gaussian field, observe 
that the expression of the 'rest' E[TM{F)f^'^\F)] appearing in p.33p is remarkably simpler 
than the ones computed in [1]. 

2. For a fixed M, the expansion (j3.33p may hold under weaker assumptions on F and /. For 
instance, if F G D^'^ has an absolutely continuous law, then, for every Lipschitz continuous 
function f :R^R, 

E[Ff{F)] = E[F]E[nF)] + E[f' (F)r^{F)]. (3.38) 
Equation (|3.38p is the starting point of the analysis developed in [9]. 
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4 Some technical estimates 



This section contains several estimates that are needed in the proof of Theorem 11.91 and Theorem 

4.1 Inequalities for kernels 

For every integer M ^ 1, we write [M] = {1, M}. Fix a set Z, as well as a vector z = 
{zi, zm) G and a nonempty set b C [M]: we denote by Zf, the element of Z^^^ (where \b\ is 
the cardinality of b) obtained by deleting from z the entries whose index is not contained in b. For 
instance, if M = 5 and b = {1, 3, 5}, then Zf, = (zi, z^, z^). Now consider the following setting: 

(a) {Z, Z) is a measurable space, and is a measure on it; 

B,q ^ 2 are integers, and bi, ...,bq are nonempty subsets of [B] such that Uj6j = [B], and 
each k £ [B] appears in exactly two of the 6j's (this implies in particular that \bi\ = 2B, 
and also that, if g = 2, then necessarily 6i = 62 = [B])', 

(7) Fi,...,Fg are functions such that Fi E L2(zI''»I,zI^'I,/xI'''I) = L2(/il'''l) for every i = l,...,q 
(in particular, each Fi is a function of \bi\ variables). 

The following generalization of the Cauchy-Schwarz inequality is crucial in this paper. 

Lemma 4.1 (Generalized Cauchy-Schwarz Inequality) Under assumptions (a) -{13) -{'y), the 
following inequality holds: 



r 1 1 
I n \PiM\l^{'izi) ■ ■ ■ l^{dzB) ^ JJ ll^i|lL2(^li>,| 

■^■^^ i=l 1=1 



(4.39) 



Proof. The case q = 2 \s just the Cauchy-Schwarz inequality, and the general result is obtained 
by recursion on q. The argument goes as follows: call A the left-hand side of (|4.39p . and assume 
that the desired estimate is true for (? — 1. Applying the Cauchy-Schwarz inequality, we deduce 
that (with obvious notation) 



A ^ ll^lllL^C^t^il) X (^i^^i a>(z,jV"^'(dZb:; 



1/2 



where the quantity 'J>(zfej) is obtained by integrating the product ni=2 ^^^"^ those variables 

Zj such that j 61. More explicitly, writing J for the class of those i € {2, g} such that bi ^bi, 

$(zbj = nii^.K)| X / n I^^Kn;,!, z,^nfef)lMl^l-l''^l(dz,c), (4.40) 

where 6^ and indicate, respectively, the complement of bi (in [i?]) and the complement of J 
(in {2, ...,g}), and /Ilje© ~ convention. By construction, one has that the sets hi such that 
2 € J are disjoint, and also that biObj =0, for every i £ J and j G J"^. If J*^ = 0, there is nothing 
to prove. If J'^ 7^ 0, one has to observe that the blocks b[ = biObl, i € J^, verify assumption (/3) 
with respect to the set [i?]\6i (that is, the class {6^ : i G J'^} is composed of nonempty subsets 
of such that Uib[ = [B]\bi, and each k G [B]\bi appears in exactly two of the b'-^s). Since 

\J'^\ ^ q — 1, the recurrence assumption can be applied to the integral on the right-hand side of 
()4.40p . thus yielding the desired conclusion. 

14 



□ 



Let be a real separable Hilbert space. The next estimates also play a pivotal role in our 
arguments. 

Lemma 4.2 Let p,q ^ 1 be two integers, r G {0, . . . ,p A q}, and u E i^®^, v S i^®*^ Then 

\\U ®rV\\fy<8>(p+q~2r) < \\u(^rV\\^'g>(p+q-2r) (4-41) 

\\u®rV\\f^®(p+q-2r) ^ ||ti||^»P Y^'ll'" ^ ll^lU^plblU®"?- (4.42) 

Moreover, if q^-\\v\\'^^^q = 1 (that is, if E[Iq{y)'^] = Ij; then 



max ||f 0^7;||^»2,-2. ^ (4.43) 



Proof. The proof of ()4.4ip is evident, by using the very definition of a symmetrized function. To 
show the first inequality in (|4.42p . apply first Fubini to get that \\u ®r ^^lly-|®(p+g-2r) = CSp-r 
u,v ®q-r v) s^S){2r) , and then Cauchy-Schwarz to get the desired conclusion. The second inequality 
in ()4.42p is an immediate consequence of Cauchy-Schwarz. Finally, the proof of (|4.43p is obtained 
by using [171 first equality on p. 183]. ^ 



4.2 Inequalities for cumulants and related quantities 

The following proposition contains all the estimates that are necessary for proving the main results 
in the paper. For every random variable Y such that i^lyl'" < oo (m ^ 1), we denote by Km^X) 
the mth cumulant of y - see Definition 13.31 Given a vector z = (zi,...,^^;) and a permutation 
a of [d], we write cj(z) = (z^^^), Zo-(^)). Given a function F{zi, Zd) of d variables and a 
permutation a of [d], we write 



(F).(z) = F(a(z)) = F(z, 



cr(l)l •••) ^o-(d) 



)• 



Also, for vectors z = (zi, Zj) and y = (yi, ...,yk), we shall write z Vy for the vector of 
dimension j + k obtained by juxtaposing z and y, that is, z V y = (zi, Zj,yi, ...,yk)- Finally, 
in the subsequent proofs we will identify vectors of dimension zero with the empty set: if z has 
dimension zero, then integration with respect to z is removed by convention. 



Proposition 4.3 We use the notation introduced in Definitions \3.3\ and \3.5i For each integer 
q ^ 2 there exists positive constants C2 (<?) , C3 (g) , C4 (g) (only depending on q) such that, for all 
F = Ig{f) with f e il®" and E[F'^] = 1, we have 



E 



^2{F) - ]^K^{F) ^ 

E\mF)w ^ 

E[\T,{F)\] ^ 



C2{q) X K4(F) 

csiq) X Ki{F) 
c^iq) X K4{F) 



(4.44) 

(4.45) 
(4.46) 
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Proof. By (j3.28p . we have slE{Ts) = Ks+i(-^) for every s ^ 1. Moreover, when F = Iq{f) is as in 
the statement, recall the following explicit representation (|3.32p : 

q [sq-2ri-...-2rs-i]Aq 

= Z]--- Yl Cg(n>---'^^)l{ri<c,}---l{n+...+r,_i<f} (4.47) 

ri=l rs=l 

X-f(.+l)g-2ri-...-2r. ( (• ■ • (/^n /)«>r2 /) • • • f)^r.j)- 

Without loss of generality, throughout the proof we shall assume that = L'^{Z, Z, where Z 
is a Polish space, Z is the associated Borel cj-field, and /x is a cr-finite measure. 

Proof of iji4-44\ )- According to (|3.28p . one has that -E'[r2(-F)] = ^K-i{F), so that the random variable 
i?[r2(-F)] — ^K3{F) is obtained by restricting the sum in ()4.47p (in the case s = 2) to the terms 
such that ri + r2 < By virtue of (|4.43p . the inequality (|4.44p will follow once it is shown that, 
for any choice of integers ri,r2 verifying such a constraint. 



||((/«'ri/)«>r2/)||^,«(35-2.i-2.2) ^ ^ max ^ 11/ (g)r /||^»2,-2.- (4.48) 

Let us first assume that r2 < q. Then ri and q — r2 both belong to {1, . . . ,q — 1}. Thus, using 
the two inequalities ()4.4ip and ()4.42p . we infer that 

I|((/®ri/)®r2/)lli3®(3?-2ri-2r2) 
«S \l\\f®q-r2 /IU»(2'-2) ll/®ri /||^0(2,-2ri) ^ _^^rnaJC_ J/ (g)^ f\\l^(2,-2r)- 

Let us now consider the case when r2 = q and ri < |. The expression 

{f®rj)®qf = {{f®rJ)J)f,m 

defines a function of g — 2ri variables. Taking into account the symmetry of / and the symmetriza- 
tion of contractions, such a function can be written as a convex linear combination of functions 
of the type 



F(t) = j /(xi,ti,w)/(x2,t2,w)/(xi,X2)d/i''+^l(w,Xi,X2) 



where w has length ri, and ti V t2 = 0"(t) for some permutation a and with t = (ti, . . . ,tg_2ri). 
Without loss of generality we can assume that ti has positive length (recall that ri < q/2 so that 
q — 2ri > 0). We denote by Sj the length of the vector xj. We then have 1 ^ si < q — ri and 
T\ < S2 ^ q — ^- Exchanging the order of integrations, we can write 

F{t) = j f{^iM,^){f®sJ) (xi,t2,w)d/i'-i+^Hw,xi). 

Squaring F and integrating, one sees that 

3 

l|-^ll|®(<7-2ri) = / n(-^®^« /)<7,(z6jd/^^(Zbl,Zb2.Zb3), 
i=l 
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with two of the Tj's equal to S2 and one to q — ri — si, where B = q + 2si, ai, i = 1,2, 3, 4, is a 
permutation of [2q — 2rj], and the sets 61,62)^3 verify property (/3), as defined at the beginning 
of the present section. It follows from Lemma 14.11 that 

3 

||-F|L®(<7-2ri) ^ max 11/ (S)r f\\l02q-2r, 

from which we deduce (|4.48|) . 

Proof of lli4-45\ )- Our aim is to prove that for any choice of (^i, r2, rs) appearing in the sum (j4.47p 
in the case s = 3 one has the inequality 

l|((/«'ri/)(ir2/)®r3/llfi8'(4<!-2ri-2r2-2r-3) ^ maX ||/ (g)r f\\j(S2q-2r. (4.49) 

l^r^tjr— 1 ^ 

Remark that {{f^rif)^r2f) already been considered when looking at T2{F) — ^k^{F), because 
of the assumption that ri + r2 < ^ . 

So, using the previous estimates and ()4.42p . we conclude directly for < q. It remains to 
consider the case when r-^ = q. 

As before, taking into account the symmetry of / and the symmetrization of contractions, it 
is sufficient to consider functions of 2{q — ri — variables of the type 

F{t) = ...,t2(g-ri-r2)) 

= / /(xi,ai,ti, w)/(x2,a2,t2, w) X 

/ (ai , a2 , t3 , X3 ) / (xi , X2 , X3 ) +'■2 (dxi , (ix2 , dx-s , dw, dai , (ia2 ) , 

where w has length ri, ai V a2 has length r2 (with either ai or a2 possibly equal to the empty 
set), and ti V t2 V = cr(t) for some permutation a. Squaring F and integrating, we claim that 
there exist integers si, S2, S3, S4 G {1, q — 1} such that 

\\Ff^2(<,-r^-r2) = j X\{f ® Si /) a, (Zfc J ^ (d^l ) • • • ^(dZfi), 
•^^^ i=l 

where B = Aq — 2(si + S2 + •Ss + •S4), Cj, z = 1, 2, 3, 4, is a permutation of \^q — 2si\, and the sets 
bi,b2,bs,b4 verify property (/3), as defined at the beginning of the present section. We have to 
consider separately two cases. 

(a): the length of X3 is not 0: we can then take si = S2 = ri and S3 = S4 equal to the length of 

X3- 

(6): the length of X3 is 0. Then either ai or a2 is not empty. Assuming that ai is not empty, we 
can take for si = S2 the length of ai and for S3 = S4 the length of X2, which is not 0. 
As before, it follows from Lemma |4. II that 

||-F||fi2{9-ri-r2) ^ JUaX ||/(g)r /||fi®25-2r, 

from which we deduce ()4.49p . 

Proof of 1^4. 4()] ). Our aim is to prove that, for any choice of (^i, r2, rs, r4) which is present in the 
sum ()4.47p (in the case s = 4), we have 

^ ^ ^ ^ 5 

ll(((/«'ri/)®r2/)^r3/)®r4/IU®(59-2ri-2r2-2.3-2r-4) ^ maX \\f ®r /||fi®2,-2r- (4.50) 

l$r$(7— 1 
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To do so, using the previous estimate (|4.45p and (j4.42p we conclude directly for < q. Hence, 
once again it remains to consider the case when = q. 

As before, taking into account the symmetry of / and the symmetrization of contractions, 
one has that the function {{if0rif)'^r2f)^r-if)^r4f is a linear combination (with coefficients not 
depending on /) of functions in 3q — 2ri — 2r2 — 2r^ variables having the form 

F(t) = t3q_2ri-2r2-2r3) 




/(xi,ai,bi,ti, w)/(x2,a2,b2,t2, w)/(x3,ai,a2,b3,t3) x 



/(bi, b2, b3, t4, X4)/(xi, X2, X3, X4)/i''+''i+'"2+''3((ix, dw, da, db), 

where w has length ri, a = ai V a2 has length r2 (with either ai or a2 possibly equal to the 
empty set), b = bi V b2 V b3 has length r^, (with some of the bj's possibly equal to the empty 
set), X = xi V X2 V X3 V X4 and ti V t2 V 13 V t4 = cr(t) for some permutation a. Squaring F and 
integrating, we claim that there exist integers si, S2, S3, S4, ,85 S {1, q — 1} such that 

5 

||F||^-,3,-2r-i-2r2-2r3 = / \\{f ® f) a,,{Zb,) KdZl) ' ' ' li{dZB), 
•'^^ i = \ 

where B = 5q — 2(si + S2 + -Ss + S4 + S5), cxj, i = 1,2,3,4,5, is a permutation of [2q — 2si], and 
the sets 61, 62) ^3) ^4, ^5 verify property (/5), as defined at the beginning of this section. We have 
to consider separately different cases. 

(a) : the length of X4 is not 0. We can then consider separately the three first factors, for which 
the same expressions as in the proof of (|4.44p are available, and the two last ones, which give rise 
to S4 = S5 equal to the length of X4. 

(b) : the length of X4 is and the length of is not 0. Then we consider separately the four factors 
which are distinct from the fourth one and proceed as in the proof of ()4.45p for them, while the 
fourth one gives rise to / (8),- /, with r equal to the length of t4. 

(c) : the lengths of X4 and t4 are 0, but the length of X3 is not 0. We then separate the five factors 
into two groups, one with /(X3, ai, a2, h^, t^) and /(xi, X2, X3), the other one with the three other 
factors. The first group gives rise to factors / <Sit /, with r equal to the length of X3, while the 
second group can be treated as in the proof of (j4.44p . 

(d) : the lengths of X3, X4 and t4 are 0, but the length of t3 is not 0. We then consider separately 
the factor /(ai, a2, b3, 13), which gives rise to a factor / /, with r equal to the length of t^. 
The four other factors can be treated as in the proof of (|4.45p . 

(e) : the lengths of X3, X4, t3 and t4 are 0. Remark that xi, X2 and b3 are non empty and, without 
loss of generality we can assume that a2 is non empty. As before, we can conclude by separating 
the five factors into two groups: for the first one we take the first factor and /(xi,X2) whereas, 
for the second one, we choose the three remaining factors. 

The desired conclusion (that is, ()4.50p ) follows once again from Lemma l4.ll 

□ 
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5 Proof of the main results 



5.1 Proof of Theorem 11.9 



The assumption = 1 implies that K := E[\F\] ^ 2. The proof follows then immediately 

from ([335]) and (OSD . 

5.2 Proof of Theorem [TTTI 

Since i?[-Fri] = and -^[-F^] = 1, one has that KiiFn) = E[F^] — 3 > 0. Moreover, because 
Fn N ~ ^(0, 1) by assumption and due again to the hypercontractivity of chaotic random 
variables, we have that Ki{Fn) = E[F^] — E[N^] — >• as n — )• oo. In the forthcoming proof we 
will need the following lemma. 

Lemma 5.1 There exists g,h £ U H C°° with bounded derivatives of all orders (except possibly 
the first one) such that E[fg'{N)] / 0, E[f'g"{N)] = 0, E[f'il{N)] = and E[f'fl'{N)] / 0. 

Proof. Let Hp, p ^ 1, denote the sequence of Hermite polynomials. Using the well-known formula 

oo ^ 

4>{x) -E[(I)^p\N)]H„{x), xeRa.e., 7V~AA(0,1), 

p=0 ^ 

valid for (/) G C°° whose derivatives are all square integrable, it is readily checked that, for almost 
all X G M, 

/c ^ ( l)^ ^ ( 1^*^ 

=cosx = > - — —H2a{x) and sinx = > — —H2a^i{x). 



(5.51) 



On the other hand, by applying several integration by parts, we can write, for h 

f'l{xf-^dx = / f^(x)(x? - \f-^dx 

-oo V^Vr 7-00 V^TT 

dxH2{x) / dy{h{y) - E\h{N)\)'^ 



oo 



1 /•+00 e~^^/^ 1 

-- j H3ix){hix) - E[h{N)])-^dx = --E[h{N)H,iN)]. 



(5.52) 



Similarly, we can prove that, for all h £U, 



E[f;i'{N)] = -\E[h{N)H,{N)]. (5.53) 

Now, let us consider g{x) = sinx. Using ()5.5ip and then (|5.52p - ()5.53p . we get that E[fg{N)] = 
7^ and E[fg"{N)] = 0. Moreover, g belongs to U because |<^"(x)| = |sinx| ^ 1, and has 
bounded derivatives. Similarly, consider h{x) = (\/e cosx — l + lH2{x)). Using once agam 
(I53T]) and then ([532]) - (1533]) . we get this time that E[fll{N)] = whereas E[fll'{N)] = -j^:^ + 
0, also with |/i"(x)| = ^ ~ ^ cosx) ^ 1 so that h ^U. 
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Proof of TheoremlTJll Recall that K3(F„) = E[F^] and Ki{Fn) = E[F^] - 3, and let g,h^Uhe 
as in the statement of Lemma 15.11 From Corollary 13.111 and Proposition 14. 3| we deduce that 

E[9{N)] - E\g{Fr.)] - \ E[f'^ {N)]ks{F,) 

^ ^\Ks{Fr,)\ \E[f'^{F^)]-E[f'^{N)]\ + ^\E[f';{Fn)]\K,{Fr,) + 2c,y^^ 

Set 

c=^ min !^l\E[f'^{N)]\,^\E[fi:'{N)]\Y 

As n ^ oo, we have E[f'J{Fn)] - E[f'^{N)] ^ 0, E[f'^'{Fn)] ^ E[f'^'{N)] = 0, and K^{Fn) ^ 0. 
Therefore, for n large enough we have that 

d{Fn,N) ^ \E\g{N)\-E\g{Fn)\ \ ^ 3c|k3(F„)| - | max{|At3(Fn)|, A^4(i^n)}. 

Similarly, we have 

E\h{N)\ - E\h{Fr,)\ -\E\r;:{N)\K^{F.^) 

^ \\k;{f^)\ \E{r;xFn)\\ + \ i^[/r(^n)] - E{ri:{N)\\ k,{f^) + 2c4h"'\\^ K,{Fnf/\ 

from which we deduce, again for n large enough, that 

d{Fn,N) ^ \E[h{N)] - E[h{Fn)]\ ^ 3cm(F„) - | max{|At3(F„)|, K4(F„)}. 

Finally taking the mean of the two previous upper bounds for d{Fn, N) yields 

3c c 
d{Fr„N) ^ y(|k3(F„)| +K4(F„)) --max{|K3(F„)|,K4(F„)} ^cmax{|K3(F„)|,K4(F„)}. 

The proof is concluded. ^ 



6 Application: estimates in the Breuer- Major CLT 

In this final section, we determine optimal rates of convergence associated with the well-known 
Breuer-Major CLT for Gaussian-subordinated random sequences - see |3] for the original paper, 
or |14j for a more modern reference. In order to be able to directly apply our previous results, we 
focus on sequences that can be represented as partial sums of Hermite polynomials. 
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6.1 General framework 



Consider a centered stationary Gaussian sequence {Xk)k£Z with unit variance and covariance 
function given by E[Xii.Xi] = p{k — I), k,l ^ X. Fix an integer q ^ 2, and set 



^ n—l 



Here, Hq stands for the gth Hermite polynomial defined by (|2.1ip . Let also N ~ ^(0, 1), and 
define '■= E\y^]^ n ^ 1. Finally, set 



, n-1 



= Y,Hq{X,). 

T '11. ^ ^ 



" fc=0 



Without loss of generality, we may assume that Xj. = X{hk), where X = {X{h) : h G S)} is some 
isonormal Gaussian process and {h^, hi)f) = p{k — I) for every k,l £ Z. We then have 



^ n— 1 

Fn = Iqifn), with /„ = —= /l^ 



nv,, 



For p ^ 1, we introduce the Banach space ^^(Z) of p-summable sequences equipped with the 
norm \\u\\iP = 

(Efcez Wik)\P)^^''. In what follows, we shall assume that p belongs to Under 
this assumption, the celebrated Breuer-Major CLT (see [3], as well |13| Chapter 7]) asserts that 

Fn ^(0, 1) as n — 7- oo. 
Remark 6.1 One has that 

Vn — ^ Piky > as n — )■ oo. 



It follows that, in the subsequent discussion, the role of the sequence Vn, n ^ 1, will be immaterial 
as far as rates of convergence are concerned. 

6.2 Explicit formulas for the third and fourth cumulants 

Let us compute, in terms of p, explicit expressions for the third and fourth cumulants of F„. 
According to Proposition 13.81 (and using the notation introduced in Section [3]), one has 

ri(i^n)=g5^(r-l)!r _ J hq~2r{fn®rfn). (6.54) 
r=l ^ ^ 



Since K3{Fn) = 2E[FnTi{Fn)] by (|3.30p . we deduce the following expression of K,^{Fn) in terms of 
the sequence /„, n ^ 1: 



2qq\{q/2 - l)!(g/2-i)^(/«' fn'^q/2fn) s^m for even q 
for odd q 
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Hence, when q is even (observe that {fn, fn®q/2fn) f^m = {fn,fn ®q/2 fn)sj»i because /„ is sym- 
metric), we have 

^3(i^n) = ^|^ E P{k-l)h{k-j)h{l-j)', (6.56) 

with d3{q) = 2qq\{q/2 — l)!(q/2-i)^- ^^^^ often useful to transform the previous expression 
as follows: we have 

Vn riy/n k,i=-j 

- -^r^ H ^"(^' - ly^'piky^'piiy^', (6.57) 

where 



/, ,s / ma,x(k,l)+ min(/c,/)_ , 
Vn{k, l)=[l ^ + ^ ) l{|fc|<„,|,|<„}. (6.58) 



Remarks 6.2 1. When q is even, one has K3(-F„) > for all n; indeed, 

n(h - l\2n(k - i\2n(l - I'^f = 

(9/2)! 



1 ^ 

E p{k-l)ip{k-j)ip{l-jY^ = —-E 



n—l /n — 1 
j=0 \fc=0 



2' 



2 



2. When q = 2, one can even prove that T2{Fn) > for all n (recall that K3(F„) = 2 £'[r2(-Fn)])- 
Now, let have a look at the fourth cumulant. Recall from |17] that 
9-1 / \ 2 



(6.59) 



For the non-symmetrized contractions of (|6.59p . the link with p is easily obtained; indeed, for any 
r = 1, . . . , g' — 1, we have: 



^ n— 1 

||/n®r-/n||J«2(,-.) = ^ E ' " •?Tp(A; " i)"-'' p{l " jT"''- (6-60) 

i,j,k,l=0 

On the other hand, we will actually face no problem due to symmetrized contractions. Indeed, 
we may forget them when deriving lower bounds, whereas we can use the inequality 

II /n®r/n ||^(82(q-r) ^ ll/n /n||^®2{9-r) (6.61) 

when dealing with upper bounds. 
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6.3 Estimates for the third and fourth cumulants 

We start with the foUowing result about the asymptotic behavior of the third cumulant of Fn- Of 
course, by virtue of (|6.55|) . only the case where q is even must be considered. 

Proposition 6.3 Assume that q ^ 2 is even. Then, 

Moreover, ifp^£^{Z), then 

^/nK3{Fn) ^ —y^V27r ^ as CO, (6.62) 

where g^^^i^) ■= p{ky^'^e^^^ is almost everywhere positive on the torus T = R\(27rZ). 

Proof. Recall the identity (|6.57p and, for any /c € Z and any n ^ 1, set Pn{k) = /9(A;)l||fc|<„} and 
\pn\{k) = \pn{k)\. Since ^ r]n{k,l) ^ 1 for all n ^ 1 and all k,l ^Z, we deduce that 



At this stage, let us recall the (well-known) Young inequality: if s,p,p' G [l,oo] are such that 
i + A = 1 + -, then 

^ ll^illfp ll^^ll^p' • (6.63) 
Hence, using first Holder inequality and then Young inequality yield 



which proves the first statement of the proposition. 

Let us now further assume that p G ^~(Z). It implies that p E i'^{Z) or, equivalently, that 
p'^/^ belongs to ^^(Z). Thus, the function g^^^ is well defined in L^(T), as being the Fourier series 
with coefficients p''/^. In particular, we deduce from Bessel-Parseval equality that 



n E /'(^^^ ^ ^ / ^9/2 ^^^^^^ as n oo. 



We also have that p^/a^^g/z belongs to that is, that Yliez p'^''^*p'^'^{l)p{iy^'^ is an absolutely 

convergent series whose value is given by Jj g^^^{t)^dt (Bessel-Parseval equality). Then, using 
()6.56p - ()6.57p and dominated convergence, we get dnSSl)- Finally, p^/^ bgi^g a covariance sequence 
as well (those of the stationary sequence --j===Hqj2{Xk)), its spectral density g^^^ is positive 

almost everywhere. (See e.g. [6].) ■ 

The situation for the fourth cumulant turns out to be not so easy, except when q = 2. 
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Proposition 6.4 There exists C > such that, for all n ^ 1, 



Ifq = 2andpe ^^/^(Z), i/i 



en 

, , LgAt^dt 
nKi{Fn) 247r -^^r^iv ^ ^ as n ^ oo, (6.65) 

where gi{t) := X^/g^ p(^)e**^* is almost everywhere positive on the torus T. If q ^ 3 and p G ^^(Z), 
i/ien liminf„_j>oo ?T''^4(-Pn) > 0. 

Proof. Thanks to ()6.59p . ()6.60p and (|6.6ip . we 'only' have to estimate, for any 1 ^ r ^ q — 1, 

^ n—l 

m = - E \Pik-l)np{^-J)np{k-^)r'■\p{l-j)rr 

i,j,k,l=0 

We immediately see that 

^ n— 1 n 
k,j=0 j=—n 

Let US first assume that g = 2, so that r = 1 necessarily. In this case, 

^(1) ^ IllPnl * \Pn\\\% ^ llPnll^, (6.66) 

where we have used Young inequality (|6.63p to get the last inequality. This proves ()6.64p when 
q = 2. Assume now that q ^ 3. By symmetry, we may and will assume that r ^ q/2. Young 
inequality (|6.63p yields that 

m ^ WKl * \Pnril ^ \\Pn\\%\\ptll- (6.67) 

This shows the desired bound when r = 1. For the other values of r (if any), we can make use of 
the log-convexity of the norms. More precisely, for a, (3 such that 2r = 2(1 — a) + {q — l)a and 
q- r = 2{l- P) + {q - l)/3, recah that 

2_^\Pn\J)\ ^ WPnWp \\Pn\\i^q-i aUQ _/__^|f'nUJI ^||Pn||£2 \\Pn\\^q-i 

We then conclude that ()6.64p holds true for any g ^ 3 as well. To finish the proof of Proposition [6? 
let us compute the limit of ?i-||/n®i/n||^^®2(g-i) under the additional assumption that p £ £^/^( 
if g = 2 and p G £^(Z) if g ^ 3. We can write 

n-l 



ri||/n«>i/n||L2(,-i) = ^ Y] Pik-iy ^p{i-jy ^p{k-i)p{l-j) 

" i,j,k,l=0 

\ E Vn{3,k,l)p{l-k)p{3)p{ky~^p{l-jY-\ (6.68) 
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where ??„(i,/c,/) = (l - ^^^''^'^ + HH£(2M_^ l{|i|<n,|fc|<n,|«|<n.} is bounded by 1 and tends to 1 
as ?i — )• oo. We know from (j6.66p if q = 2 and from (j6.67p if q ^ 3 that the series 

p{i-k)p{j)p{kr^p{i-jr^ 

j,k,l&Z 

is absolutely convergent under our assumption on p. By dominated convergence, we get that 
n||/n<»i/n|||®2(9-i) ^ {p* P'^^^)'^ (0 = \\P* P''~^\\% as n ^ oo, 

implying in turn that liminf„_j.oo ?T-«^4(-?*n) > as expected. Finally, when q = 2 we have Vn — )• 
- I gi{tfdt, so that ([6:65]) holds true thanks to ([639]) . 



6.4 Breuer-Major CLT 

The following result is the so-called Breuer-Major theorem (it is called this way in honor of the 
seminal paper [3]). For sake of completeness, we provide here a modern proof, that relies on the 
bounds (j6.64p for ^^{Fn) and Theorem ll.il 

Theorem 6.5 (Breuer-Major) Assume p £ £'?(Z). Then converges to ^(0,1) in total 
variation as n —t- oo. 

Proof. By virtue of Theorem 11.11 it is (surprisingly) enough to prove that Ki{Fn) tends to as 
n —7- oo. Let us first assume that q = 2. Then, using Holder inequality for k > M, we have, for 
any n > M, 

3 / \ 2 



n-l 



\k=-n+l J \\k\^M I \J\/<|A:|<ri 



2 



We can then conclude by a standard argument (choosing M large enough). For g ^ 3, we proceed 
analogously, after having noticed that 

2 2 
n — 1 \ / n— 1 \ / ^ n — 1 \ / -i 1 



with 

^ E 1^(^)1'^' ^ ( E 1/^(^)1' 1 ^ '^"'\\p\\V\ 



\fc=— n+l / \A,'=— n+1 / \ k=—n+l / \ k=—n+l 



1 / X 1-1/9 

n— 1 ' ^ 



k=-n+l \|fc|<n 



and, for any M < n, 

E wm^*;^ E + E 

Hence, by choosing M large enough, we get that K.4,{Fn) tends to 0, and the proof is concluded. 
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6.5 The discrete- time fractional Brownian motion 

Let Bh be a fractional Brownian motion with Hurst index H G (0,1). We recall that Bh = 
{BHit)}t£R is a centered Gaussian process with continuous paths such that 

E[BH{t)BH{s)] = + -\t- s,teR. 

The process Bh is self-similar with stationary increments, and we refer the reader to Nualart |16| 
and Samorodnitsky and Taqqu |20| for its main properties. In this section, we offer fine estimates 
for K3{Fn) and K4(-F„) respectively, when the sequence {Xk)ke:Z corresponds to increments of B^ , 
that is, 

Xk = Bnik + l) - Bnik), k £ Z. (6.69) 

The Xjg's are usually called 'fractional Gaussian noise' in the literature, and are centered stationary 
Gaussian random variables with covariance 

pik) = E[XrXr+k] = ^ (1^ + M^" - 2|A:P^ + \k- l^^) , r, fc G Z. 
The covariance behaves asymptotically as 

p{k) ^ H{2H -l)\k\^^-^ as|A:|^oo. (6.70) 
In particular, when H > I we observe that, for \k\ large enough. 



2 

1 



pik)^HiH--)il + \k\r^-'. (6.71) 

In the following results, we let the notation introduced in Section [6.11 prevail, and we assume 
that the sequence {Xk)k& is given by ()6.69p . We also use the following convention for non- negative 
sequences {un) and (f„) (possibly depending on q and/or H): we write Vn ^ Un to indicate that 
< liminf.„^oo ■Vn/ii„ ^ lim sup^^^^ < oo. 



Proposition 6.6 Assume that q ^ 2 is even. We have: 

3q 



n-2 o<iJ<l-|- 



1 



n 2 log^ n if H = I — 



Proof. When i7 < 1 - ^, we have that p £ so that 

< liminf "v/n K3(F„) = limsup \/n'^3(-fn) < oo 

71— >oo n— 5>oo 

by Proposition 16.31 For the other values of consider first the limsup. We have that 



o / . I loff n if H = 1 — 

\p{k)\"^/' ^ c \ 



26 



from which we deduce the finiteness of the hmsup by Proposition 16.31 together with the fact that 
Vn q\ Ekez Pi^y ^ (0' o«) when H <1-^. 

Let us now focus on the liminf. Since > 1 - ^ > i, the lower bound ([GTT]) holds for 
\k\ large enough. In fact, by considering a decomposition of the type p = + t with r having 
compact support, we can do as if the inequality (|6.7ip were valid for the small values of as 
well. Indeed, it suffices to use the fact that Ks{Fn) comes from a trilinear form (see ()6.57p ) as well 
as the inequality ^ ^ (jxl'? + , so to deduce that, when there is at least one term r, 

the asymptotic order of the corresponding contribution in ^^{Fn) is 0(n~2). Recall the definition 
(j6.58p of r/„(fc, I). With the extra assumption that (j6.7ip holds true for all k, we get that 

Y,rin{hl)p{k-lYl^p{kfl^p{lfl^ 

k,l&Z 

is bounded by below by a constant time 

^ ^ (l + |A;-Z|)(^^-ik(l + |A:|)(^-i)'?(l + |/|){^-i)''. 

Note that, for k,l with |/| < \k\ ^ 2\l\, we have l + \k-l\ ^ 2(1 + \k\) as weh as 1 + |A; - /| ^ 
3(1 + |/|), so that 

which, by summing first over k and then over I, concludes the proof for the liminf. ■ 

By reasoning similarly, we obtain an estimate for K4(-Fn)- (In the following statement, the 
limsup is partially known from [2].) 

Proposition 6.7 For q G {2,3}, we have 
1 

-1 



KiiFn) > 

whereas, for q > 3, 
K,i{Fn) > 



AqH-Aq+2 



= 1 - ^ 



H 
1 - 



(6.72) 



-1 

1 log(n) 

AH-A 



AH-A 



log n 



AqH-Aq+2 



< F < f 
1<H<1 



H 
1 - 



1 



2g-2 



(6.73) 



2q~2 



2q-2 

<H <1 



j_ 

2q 



Proof. The proof of the limsup is straightforward by using Proposition l6.4l To get the liminf result, 
we first consider p = + r as in the proof of Proposition 16.61 and verify that the contribution 
of all the terms containing at least one r are of lower order. Here again, we can therefore do as 
if the inequality (j6.7ip were valid for all k. This allows us to bound (|6.68p by below by following 
the same line of reasoning as in the proof of Proposition 16.61 We finally conclude thanks to (|6.59p . 
Details are left to the reader. ■ 

By comparing the asymptotic behaviors of K3(-F„) and K4(F„), we observe the following non- 
expected fact. 
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Corollary 6.8 When q ^ 6, there exists a non-trivial range of values of H for which K4{Fn) 
tends less rapidly to than K^{Fn). 

Proof. The two functions that give the behavior of n-s^Fn) and K4(F„) are piecewise hnear and 
concave. It is therefore sufficient to consider the value H = 1 — For this value, one has 
f^AiFn) ^ n""-^ » n-1/2 when ^ < ^ < i. ' ■ 
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